Extreme statistics of complex random and quantum chaotic states 
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An exact analytical description of extreme intensity statistics in complex random states is de- 
rived. These states have the statistical properties of the Gaussian and Circular Unitary Ensemble 
eigenstates of random matrix theory. Although the components are correlated by the normalization 
constraint, it is still possible to derive compact formulae for all values of the dimensionality A*'. The 
maximum intensity result slowly approaches the Gumbel distribution even though the variables are 
bounded, whereas the minimum intensity result rapidly approaches the WeibuU distribution. Since 
random matrix theory is conjectured to be applicable to chaotic quantum systems, we calculate the 
extreme eigenfunction statistics for the standard map with parameters at which its classical map is 
fully chaotic. The statistical behaviors are consistent with the finite-A formulae. 
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The study of the statistics of extreme values [ij has 
found many applications in diverse areas such as geo- 
physics, metereology, economics, structural engineering, 
ocean waves, and dynamical systems. The subject is cur- 
rently undergoing a resurgence of interest due to recent 
catastrophic events such as hurricanes, floods, and a par- 
ticularly deadly tsunami as well as a number of research 
advances 0]. The kinds of questions being asked are, for 
example, what are the distributions for extreme events, 
or what are the inter-event gap distributions. It has long 
been known that if the underlying events are indepen- 
dent and identically distributed, then for appropriately 
rescaled variables there are three possible limiting uni- 
versal distributions for the extreme maximal events: the 
Frechet, Gumbel and Weibull distributions [ij. Respec- 
tively, they arise depending on whether the tail of the 
density is a power law, or faster than any power-law, and 
unbounded or bounded. If there are correlations, then it 
is known that these universal distributions are reached 
for a sufficiently fast decay of auto-correlations Q . 

In this letter, we apply these powerful methods to 
study the extreme properties of random vectors, or wave 
functions more generally, in the context of quantum me- 
chanics. Our motivation is that the eigenstate intensi- 
ties in fully chaotic systems with no particular symme- 
tries are conjectured to behave exactly as these random 
vectors subject only to a normalization constraint as in 
random matrix theory. For chaotic systems, the applica- 
bility of random matrix theory [3, Q has been well ap- 
preciated since the Bohigas-Giannoni-Schmit conjecture 
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which states that the spectral fluctuations of quantized 
classically chaotic systems can be modelled by a suit- 
able ensemble of random matrices @. In fact, a certain 
number of extreme spectral properties have already been 
derived 0, H, 0, beginning with the well-known result 
for the distribution of the largest eigenvalue [ll|, [13] ■ The 
corresponding treatment of random vectors or quantum 
eigenvectors has not yet been addressed. However, see 
[l3| for an initial foray into random waves. 

In fact, eigenstate intensities in strongly chaotic sys- 
tems are known to follow an exponential density, which is 
consistent with states uniformly distributed over a stan- 
dard simplex ij], as happens in the Unitary Ensem- 
bles (if an anti-unitary symmetry is respected, Porter- 
Thomas density, hypersphere, and Orthogonal Ensem- 
bles A similar class of problems shows up in frag- 
mentation, i.e. randomly cutting an object of fixed length 
into N pieces [l6| . 

It is possible to give compact, exact formulae for all 
dimensionality N in spite of the correlations introduced 
by the normalization constraint. It turns out that the 
small- A*" distributions for the maxima differ considerably 
from their asymptotic limit (which turns out to be Gum- 
bel, exp[— e""^]) thus giving the possibility of extracting 
system size from the distributions. In an A^-dimensional 
complex Hilbert space a general state is represented in 
a fixed orthonormal basis \i) as = X^i^^i ^iN)- If the 
Zi are complex components of a random state, they are 
not constrained by any requirement other than normal- 
ization, and their joint probability distribution is: 



P(Z1,Z2, . . . ,Zn) 
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The real and imaginary parts of the components are 
spread in an unbiased, microcanonical, manner over the 
2A^-dimensional unit sphere. The intensities jz^p are dis- 
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tributed uniformly on an iV — 1 simplex. Thus a random 
state is correlated by this requirement. Given such an 
ensemble of random states we consider the probability 
distribution p{t,N) of t = max{|zi|^, |z2|^, • • • , I^Arp}. 
Let F{t,N) be the probability that all \zj\'^ < t. This 
is called the distribution or cumulative density of such 
extreme events, i.e. F'{t,N) ~ p{t,N) where the prime 
denotes differentiation with t. Thus 
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where rj = \zj\. Use a Fourier decomposition of the 5- 
function and perform the {9i, r^} integrals to arrive at 
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Next expand the power in a binomial series to find 



(3) 



N 
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where 
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cj^-isign(cj). (5) 
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The exact result for F{t,N) is 

^6) 

Expanding the power again and using the identity [l7|: 
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gives "resummed" expressions for F(t, N) valid in the 
intervals 4 = [l/(fc+l), 1/k], where A: = 1, 2, • • • , iV - 1: 



Fiteh,N)=Y, 
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and F{t < 1/N,N) — 0. Thus the cumulative density is 
a piecewise smooth function on the intervals Ik • 

Given the simple form of the distribution above, it is 
useful to interpret them combinatorially and derive them 
from such an approach. First note that 
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Pl{zi,Z2, ...,Zl) = 
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where Pi is the reduced probability density for I complex 
components (P/v is identical to the full distribution P 



above). If < € /i = [1/2,1] there can be at most only 
one such component. Therefore, the fraction of states 
with a component larger than t is exactly the fraction of 
components larger than t. Since the desired quantity is 
the fraction of states such that all components are less 
than t, it is the simply the complement: 



F{t,N) J 
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Pi{z)dz 
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The factor N accounts for multiplicity of choice of this 
one component. This integral is elementary for the com- 
plex case and gives F{t,N) = 1 — A^(l — t)^^^, which 
agrees with the series in the RHS of Eq. ([5]) which termi- 
nates at fc = 1 for t S /i . 

If t G I2, it is possible that there are at most two 
such components. The number of components > i is no 
longer the number of sequences (states) with at least one 
component more than t as it double counts states which 
have two components larger than i, a contribution which 
must get subtracted. This same logic extends, and in 
the next interval /a, the number of pairs over counts the 
contributions of triples. Similarly, this reasoning carries 
forward to any distribution with a unit norm constraint 
and gives for the cumulative density 

F{telk.N)^ 

l^K-l)" / Pm{xi,...,Xm)dxi---dxm, 

^0 J\.,\^>t 
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which is the generalization of Eq. ([8]) (the m — term in 
the above expression is unity) . This generalization could 
be used as the starting point for an analysis of real ran- 
dom states as well as general density matrix eigenvalues 
whose sum is also constrained to be unity. It is interest- 
ing that piecewise continuous extreme distributions have 
been identified in dynamical systems [3, [l^ and frag- 
mentation problems |l6|, where a similar combinatorial 
approach was also applied. 

From the distribution Eq. ([6]) above, it is possible to 
derive exact formulae for the moments. In particular, 
results for the first (mean) and the second moments of 
the maximum components are 
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where H{N, k) is the Harmonic number of order k defined 
by the finite sum X]m=i ^ind 7 = 0.5772 • • • is the 

Euler-Mascheroni constant. The standard deviation can 
be calculated exactly and its large- form is: 
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FIG. 1: The convergence of the exact probabihty density to 
the asymptotic Gumbel distribution using the scaled variable 
X = N{t — ln{N)/N) with increasing A^. The inset shows the 
difference between the exact and the Gumbel distribution for 
the same values of N, but in the unsealed variable. 



It turns out that the asymptotic N oo hmit here is 
also the result for uncorrelated exponentially distributed 
variables of mean 1/N. This limit may be calculated by 
simply taking Pi as the probability density of an inde- 
pendent process, which gives 
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Expressed in terms of the standard linearly scaled vari- 
able X = {t — ajq) /bjq , this distribution is seen to co- 
incide with the Gumbel distribution where the param- 
eters are given by = \n{N)/N and = 1/N. As 
should happen, the mean and the standard deviation cal- 
culated from the Gumbel distribution coincide with the 
leading order contributions derived in Eqs. (|12ll4p . That 
the maximal statistics is a Gumbel distribution is inter- 
esting because the intensities have a finite support, being 
constrained to lie on the standard N — 1-simplex. The 
expectation for uncorrelated bounded variables would be 
the Weibull distribution; however the correlations intro- 
duced by the normalization constraint transforms this to 
the Gumbel distribution. In Fig. ([T]) we compare the ex- 
act probability density with that of the Gumbel density 
(exp[— a: — e~^]) after appropriately rescaling. It is clear 
that the approach of the exact density to the Gumbel one 
is rather slow and at N around 100 there are still signifi- 
cant differences between the asymptotic and the exact. It 
is also instructive to note that without the rescaling the 
densities actually diverge as N increases. This difference 
gives the possibility of extracting system size information 
from the extreme statistics. 

The distribution of the minimum intensity s = 



min{|zi| 



^2 



|zAr| } on the other hand is a much 



tribution. The fraction of states such that the minimum 
is larger than some s is the fraction of states such that 
all the components are larger than s. Thus if F{s, N) is 
the cumulative distribution of the minimum, it is given 
by 



F{s,N) = 1-- 



(16) 

which is very similar to the integral in Eq. Its eval- 
uation proceeds similarly to the maximum, and gives: 



F{s,N) = 



1 - (1 - Ns)^-'^ < s < 1/A^ 
1 1/N <s<l 



(17) 



It is clear that the minimum cannot exceed 1/A^, just as 
the maximum cannot be less than this. The average min- 
imum component is easily calculated and is exactly equal 
to (s) = I/N"^. The distribution for the minimum does 
not have the piecewise continuous character we observed 
for the maximum. This has a geometrical interpretation 
in terms of the standard-simplex. In the case of the max- 
imum the integral in Eq. ([2]) can be interpreted in terms 
of volumes of subsets contained in the region bounded 
by the standard N — 1-simplex; these volumes enclosing 
more complex shapes for increasing t, as they pierce the 
simplex boundary. Whereas for the integral in Eq. (|16p 
the volumes involved are those of the entire simplex and 
the volume of a subset that never pierces the A'^ — 1 sim- 
plex. 

For large N the distribution of the minimum ap- 
proaches the exponential: 



F{s,N oo)^l- exp(- 



-N\s) 



(18) 



simpler quantity and is not asymptotically a Gumbel dis- 



This being a special case of the Weibull distribution, it is 
indeed what one would expect of uncorrelated variables 
with a compact support. That the minimum cannot be 
less than zero presents a strong constraint and for small 
components the normalization correlation is not so im- 
portant. Thus, the distribution of the maximum and 
minimum of the complex random vector intensities fol- 
low different universal distributions asymptotically. It is 
noteworthy, however, that the limiting large deviations 
of the maximum component toward a small value which 
occur in the interval In-i = ^/N < t < 1/{N — 1) is 
distributed as F{t, N) — {Nt — 1)^^^ which is an ex- 
act reflection about the value 1 /N of the behavior of the 
minimum component. 

In order to compare these extreme statistics to 
the statistical properties of the eigenfunctions of a 
Hamiltonian system, consider a quantum kicked rotor 
[20| . This is a stroboscopic mapping of a kicked one- 
dimensional particle of unit mass moving on a circle of 
unit perimeter with the Hamiltonian H{q,p) = p'^/2 — 
(if/4^2) ^og(2^g) 5(t - n). From Hiq,p), the 
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x=Nt-ln N 

FIG. 2: The probability densities (histograms) of the scaled 
maximum and minimum (inset) intensity of eigenfunctions in 
the position basis of the quantum kicked rotor for A'^ = 32 in 
the parameter range 13.8 < K < 14.8. Shown as a continuous 
line is the exact density for random states while the dotted 
ones are the respective Gumbel and WeibuU densities. 

classical mapping can be given 21] (modulo unity): 

Pi+i =Pi- K/(2it) s\n{2nqi) ; g^+i = + Pi+i (19) 

where the potential kick is applied before the free motion 
and the phase space is restricted to the unit torus {q,p £ 
[0, 1)) for convenience. For iiT ^ 5 the map is highly 
chaotic, although the phase space almost always contains 
some tiny proportion of regular motion mixed in. 
In the position basis, the evolution operator is 

1 fiNK /2n(n' + a]\\ 
^""' = iV^"H^'°'l N )) 

m=0 ^ ' 

(20) 

The two phases < a,/? < 1 can be used for control- 
ling parity and time-reversal symmetry breaking (a = 
1/2, /3 = preserves both symmetries). Choosing /? well 
away from or 1/2 breaks time-reversal invariance, which 
for if ^ 5 leads to quantum chaotic states that are 
complex and whose extremes should follow the distribu- 
tions of the unitary ensembles. The dimensionality of the 
Hilbert space N is the inverse (scaled) Planck constant. 

An ensemble of roughly 30, 000 quantum chaotic eigen- 
states is created from the N orthonormal states, as well 
as from those obtained by a variation of parameters if, 
a, and (5 such that while the quantum spectrum is sig- 
nificantly changed, the classical dynamics is not. Fig- 
ure ^ shows the density of the maximal and minimal 
intensities in the position basis for a range of K values 
where the classical map has no significant islands and 



is highly chaotic [22|. The derivatives of the exact re- 
sults in Eqs. (|8ll7p fit the quantum system histograms 
very well. The deviations are roughly of the scale of the 
expected sample size errors. Note that the dynamical 
system results require the exact density for the maxi- 
mum as the asymptotic approach is too slow. This is in 
contrast to the exact finite- TV Dyson/Mehta fluctuation 
measures [3], which approach their asymptotic limits so 
quickly that finite- results are rarely, if ever, used. The 
inset shows that the fit to the WeibuU is excellent even 
for the small value of TV used. Given the excellent agree- 
ment of the dynamical systems results with the analytic 
forms, deviations may in turn be used to investigate the 
important issue of eigenstate localization. 

In summary, recent years have seen a fast and fruitful 
development of the merging of extreme statistics tools 
and random matrix eigenvalues. This letter extends this 
general perspective to the properties of eigenvectors. 
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